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Abstract 

(~| , Scalar field dynamics may give rise to a nonzero cosmological variation of 

^ i' fundamental constants. Within different scenarios based on the unification 

of gauge couplings, the various claimed observations and bounds may be 
' combined in order to trace or restrict the time history of the couplings and 

, masses. If the scalar field is responsible for a dynamical dark energy or 

quintessence, cosmological information becomes available for its time evo- 
, lution. Combining this information with the time variation of couplings, 

^ ' one can determine the interaction strength between the scalar and atoms, 

OO , which may be observed by tests of the Weak Equivalence Principle. We 

(N . 

compute bounds on the present rate of coupling variation from experi- 
, ments testing the differential accelerations for bodies with equal mass and 

• ' different composition and compare the sensitivity of various methods. In 

, particular, we discuss two specific models of scalar evolution: crossover 

QQ ' quintessence and growing neutrino models. 

o ■ 

■ 1 Introduction 

X. 

^ , Several positive and null results on the time variation of fundamental constants 

have been reported. In a previous paper [Ij we reviewed the relevant obser- 
vational evidence and bounds, and introduced several distinct scenarios with 
unification of gauge couplings (GUT) [2]. Each scenario is characterized by 
specific relations between the variations of different Standard Model (SM) cou- 
plings: specifically, all fractional variations in SM couplings were taken to be 
proportional to the variation of one underlying unified coupling or a fixed linear 
combination of several unified couplings. The different scenarios have different 
consequences for the interpretation of any signal of nonzero variation. We also 
divided the observations into several cosmological "epochs" and discussed the 
consistency between them within each epoch, as well as whether a monotonic 
time evolution between epochs could fit the observations. Up to this point no 
particular dynamical mechanism for the coupling variations was considered. 

In quantum field theory, a possible time variation of couplings must be as- 
sociated to the time variation of a field. This field may describe a new "funda- 
mental particle" or stand for the expectation value of some composite operator. 
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In this paper we discuss the consequences of the simplest hypothesis, namely 
that the field is a scalar, such that its time varying expectation value preserves 
rotation and translation symmetry locally, as well as all gauge symmetries of 
the Standard Model. 

Any observation of coupling variations in the cosmological history from nu- 
cleosynthesis up to now would imply that the expectation value of this scalar 
field changes appreciably during this cosmological epoch. Such a "late time evo- 
lution" (in contrast to the electroweak or QCD phase transition, GUT-phase 
transition or even inflation) is characteristic for models with a dynamical dark 
energy or quintessence. Indeed, the potential and kinetic energy of the scalar 
"cosmon" field would contribute a homogeneous and isotropic energy density 
in the Universe, and therefore lead to dynamical dark energy [3j. The time 
variation of "fundamental couplings" is a generic prediction of such models [3]. 
In fact, it has to be explained why the relative variation is tiny during a Hubble 
time, while the simplest models for a scalar cosmon field arising from super- 
strings or other unified frameworks typically lead to relative variations of the or- 
der one. Efficient mechanisms have been proposed for explaining the smallness 
of variations, for example based on the approach to a fixed point [5l[6]. Unfor- 
tunately, the strength of such mechanisms is not known, such that the overall 
strength of the scalar coupling to matter remains theoretically undetermined. 
Several investigations have explored the consequence of scalar quintessence 
models for the time variation of couplings, as well as models specifically designed 
to account for coupling variations [3[3[9l[ini[IIl[I3[Il[Il[l5l[l6l[l8l[I7l[l9]. 

The cosmon coupling to atoms determines both the outcome of tests of the 
Weak Equivalence Principle and the time variation of couplings in the most 
recent cosmological epoch, including the present. For the latter, one also needs 
information for the rate of change of the expectation value of the scalar field, 
which can be expressed in terms of cosmological observables, namely the fraction 
in dark energy, Qh, and its equation of state, Wh- In consequence, the differen- 
tial acceleration r] between two bodies of different composition, and therefore 
different specific "cosmon charge" , can be related to the present time variation 
of couplings and cosmological parameters as [13j 



Here the present dark energy fraction is ~ 0.73 and ct/a is the present 
relative variation of the fine structure constant. The "unification factor" F 
encodes the dependence on the precise unification scenario (defined in Section 
[2]) and on the composition of the test bodies. In the notation of [13] we have 
F = ARz{l + Q), where ARz = A{Z/{Z + N)) is the difference in the proton 
fraction Z/A between the two bodies, which takes a value ARz < 0.05 for 
typical experimental tests of the equivalence principle. The present dark energy 
equation of state Wh is close to —1; we take 1 + Wh ^ 0.1. The factor F will 
be calculated for our different unified scenarios in Section 14.21 for typical test 
mass compositions we find 1 < < few x 10^. 

Once F is fixed, the relation ([T|) allows for a direct comparison between 
the sensitivity of measurements of rj versus the measurements of a /a from 
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laboratory experiments, or bounds from recent cosmological history, such as 
from the Oklo natural reactor or the isotopic composition of meteorites. 

Our paper is organized as follows. After a brief summary of our previous 
results concerning evolution factors for different unified scenarios, we introduce 
in Section O a cosmologically evolving scalar field as the source of the dynamics 
of a possible coupling variation. The scalar has both an effect on large-scale cos- 
mological evolution as a consequence of its homogeneous potential and kinetic 
energy, and may also produce local gravitational effects, violating the Weak 
Equivalence Principle (WEP) due to its interactions with matter. In Sec. [4] 
we use both aspects to set further constraints on a possible time variation of 
couplings. We determine, for each unified scenario, which observational meth- 
ods are most sensitive to detect late-time or present-day variations. This may 
enable us to distinguish between scenarios if there is a non-zero signal. 

In Section [5] we determine for each scenario the evolution factors that a 
smooth time evolution of the scalar must satisfy to be consistent with observa- 
tions, in particular with a detectable nonzero variation at z ~ 1-3. We discuss 
whether a simple model of monotonically varying "crossover quintessence" can 
produce the required behaviour. Section [6] is devoted to a recently proposed 
type of model [201 EI] where the scalar field evolution is halted by its coupling to 
neutrinos whose mass increases. Such models solve the "coincidence problem" 
of dark energy and give rise to an interesting scalar field evolution with damped 
oscillations at recent epochs. These models need not obey bounds derived for a 
monotonic time variation, and we use instead global fits over all data for some 
specific choices of cosmological model. In Section [7] we summarize and discuss 
our findings. 



2 Unified scenarios and epochs 

In [1] we defined different GUT scenarios for varying couplings. They are based 
on the unification of gauge couplings at a scale Mx below the Planck mass 
Mp, with a unified coupling ax- We considered also variations of the Higgs 
v.e.v. {(j)) and, for SUSY models, of the superpartner mass scale m. Under the 
hypothesis of proportionality of relative variations of the "unified couplings" 
ax, Mx/Mp, ((/))/Mx, fh/Mx we can write 

Aln^ = dM/, ^\nax = dxl, Aln^ = dj,/, /\\n^ = dsl. (2) 
Mp Mx Mx 

The "evolution factor" l{z) depends on the redshift z and will be linked to the 
evolving value of a scalar field. If ax varies nontrivially we may normalise I 
via dx = 1. In models without supersymmetry we set ax = 1/40 and dg = 0, 
while for supersymmetric theories we take ax = 1/24. As explained in [Ij, we 
neglect variations of Yukawa couplings, thus all masses of the charged fermions 
in the Standard Model vary proportional to (0). 

The different unified scenarios investigated in p!] can be characterized by 
the unification coefficients d^ as follows: 

• Varying a alone 
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• Scenario 1: Varying gravitational coupling, dM = ^, dx = dn = ds = 

• Scenario 2: Varying unified coupling, dx = 1, dM = dn = ds = 0: we 
distinguish Scenario 2 (non-SUSY) and 2S (SUSY) if relevant 

• Scenario 3: Varying Fermi scale, dn = 1, dpj = dx = ds = 

• Scenario 4: Varying Fermi scale and SUSY-breaking scale, dff = ds = 1, 
dM = dx = 

• Scenario 5: Varying unified coupling and Fermi scale, dx = 1, dn = 7, 
dM = ds = 0: we considered 7 = 42 for illustration 

• Scenario 6: Varying unified coupling and Fermi scale, dx = I, dn = ds = 
7, dM = 0: we considered 7 = 70 and 7 = 25. 

Computing the couplings of the Standard Model of particle physics in terms of 
the unified couplings, all relevant observables are then found in terms of I, dM, 
dx, dn and ds in each scenario. This allowed us to reduce each observational 
result to a constraint on the variation of the underlying unified parameter, or 
equivalently on l{z). 

For each unified scenario, there are a number of questions that may be 
addressed within our framework. First, how many of the claimed signals of non- 
zero variation may be consistently accomodated, given a monotonic variation 
and given other observational constraints on variations. 

As already stated, atomic clock experiments and WEP violation put bounds 
on the rate of time variation at present. These bounds can be extrapolated to 
the past and compared to observational constraints at previous epochs, given 
some assumption for the form of time variation. For this comparison we will 
generally assume a constant time derivative and restrict the extrapolation to 
small redshifts. One may also consider the case of oscillating (non-monotonic) 
variation. However, this introduces too much freedom to allow a general anal- 
ysis. Hence we consider a few specific models with oscillations in Sec. [H 

Second, we may consider specifically BBN. If a nonzero variation occurs at 
high redshift, as suggested by spectroscopic determinations of a and/or /i, one 
can ask whether such variation can improve the fit to observationally determined 
abundances including ''Li. If a nonzero variation at BBN is favoured, what is 
the range of values for the evolution factor Iq = Ibbn, and what is the best fit 
that can be obtained for the BBN abundances? 

Third, what additional signals are expected for a moderate increase in ob- 
servational sensitivity, particularly in the case of tests that at present give null 
results? In other words, what method could rule out or confirm any given 
model most efficiently? In section 14.41 we consider, in particular, probing cur- 
rent claimed variations by atomic clocks or detecting the corresponding scalar 
couplings by WEP experiments. 

In order to simplify the comparison of different observations we defined 
several "epochs" as follows: 

• Epoch 1: Today until Oklo, z < 0.14 
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• Epoch 2: Cosmological probes at low redshift, 0.2 < z < 0.8 

• Epoch 3: Intermediate redshift, 0.8 < z < 2.4 

• Epoch 4: High redshift, 2.4 < z < 10 

• Epoch 5: CMB, z « 1100 

• Epoch 6: BBN, z ^ 10^°. 

Laboratory bounds on the present-day time variation of fundamental couplings 
are incorporated into Epoch 1 by simple linear extrapolation back to the time of 
Oklo, 1.8 billion years ago. Then for each scenario we evaluated a weighted aver- 
age of the observations in each epoch, resulting in "evolution factors" li, . . . , Iq. 
The relevant values of evolution factors k are listed in Table 4 of [T]. For con- 
venience, we reproduce this table in the Appendix of this paper. 

We found that no scenario was able to fit all the nonzero signals or hints 
of variation that we considered with a monotonic time variation. The possible 
nonzero variations are the Murphy et al. observations of a variation in 143 
quasar absorption systems [22j; the Reinhold et al. observations of fi = nip/me 
variation from H2 spectra at high redshift [23], and the potential solution of the 
^Li problem of primordial nucleosynthesis (BBN) via variation of quark masses 
[241 1251 [26] . However, we were able to obtain monotonic variation within some 
scenarios, if we considered that one of these nonzero results was affected by an 
unknown systematic error that increased the uncertainty. For example Scenario 
6 with 7 = 25 is consistent with monotonic evolution if the more conservative 
yu bound of [27] is substituted for the value of /i variation obtained in [23] . 

3 Coupling variations from evolution of a 
cosmological scalar field 

We now consider the dynamics of the degree of freedom that gives rise to the 
variations of couplings. If we simply wrote down an action with the couplings 
ax, ete., promoted to functions of cosmological time, then almost all the sym- 
metries that have been used to determine the form of underlying theories would 
be explicitly violated. This would leave well-tested properties like Lorentz in- 
variance, charge conservation and energy-momentum conservation without a co- 
herent theoretical explanation, and would destroy any predictivity since physics 
would depend on arbitrary, unknown functions of time. We must introduce one 
or more new propagating degrees of freedom, that allow us to take a limit where 
the fundamental symmetries are exactly restored if the corresponding fields have 
static values. This will lead to a closed system of dynamical equations giving, 
in principle, predictions for the cosmological and local/present-day behaviour 
of couplings. 

The simplest possibility of a new degree of freedom allowing a limit of local 
Lorentz invariance is a scalar field, which may be fundamental or effective (com- 
posite). If the scalar is a gauge singlet, all gauge symmetries and the associated 
charges are conserved. In cosmology we may consider a homogeneous field that 
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only depends on time. More precisely, the value of the field is statistically ho- 
mogeneous on large scales over hypersurfaces on which the matter density is 
also statistically homogeneous. For a very light field weakly coupled to matter 
the local perturbations are generally small relative to the cosmological evolu- 
tion. In other words, the evolution of the scalar field in a cluster of galaxies 
or on Earth does not decouple from the cosmological evolution (in contrast to 
the gravitational field), such that its cosmological time evolution is reflected in 
a universal variation of couplings, both on Earth and in the distant Universe 
[m [281 [29]. 

Observable coupling variations are proportional to the interaction strength 
of the scalar field to the constituents of ordinary matter (atoms). In the limits 
of either weak interaction or a slow time evolution of the scalar the coupling 
variations may become tiny. It is tempting to identify the scalar with the 
cosmon that is the key ingredient for dynamical dark energy or quintessence, 
but we will first keep the discussion general. In any case cosmology will put 
additional constraints on the evolution of such a field. 



3.1 Coupling functions 

The assumption of unification implies that the variation of the unified coupling 
is related via some coupling function Bx to the scalar ip 

ax = Bx{(p), (3) 

realized typically by a gauge-invariant kinetic term for the gauge fields in the 
GUT-model oc Bx{'^)~^F^^'^ F^y [311[32].0 In a linear approximation, this yields 

B' 

Alnax = Aln5x('/') - TT^^V^ = c^x^ (4) 

Bx 

where the dimensionless scalar ip is obtained by dividing out the reduced Planck 
mass Mp = (SvrG)"^/^. Similarly, the unification mass Mx{<p) or the Fermi 
scale {(p){(p) may depend on y?, with 

Aln^ = AlnBM{ip), A In = A In ((/,). (5) 

More generally, one has in the linear approximation for the couplings Gk of the 
Standard Model of particle physics 

Aln Gk = dkl ^ Pk^V, (6) 

where we define 

(5k = -^lnBk. (7) 

For Px 7^ we may choose dx = 1 and I = (3xAip. (As mentioned in Section 
[21 in scenarios where ax does not vary we take dn = 1, I = PnAip.) All 



^For the pure electromagnetic field the gauge invariant description of a varying fine struc- 
ture constant was given in [331 134j . 
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quantities may depend on redshift or the epoch n. For the relation ([2]) to hold 
it is sufficient that the ratios AlnB^/AliiBx = dk/dx remain constant. 

In this and the next two sections we investigate the hypothesis of monotonic 
evolution. This assumes first a monotonic evolution of the scalar field towards 
today's value ipQ. The presently measured laboratory values of the fundamental 
'constants' correspond to = ipQ. Secondly, the coupling functions Bk{(p) are 
also assumed to be smooth and monotonic. Thus 

ip{z) -^Q = Av3(z) (8) 

and l{z) are monotonic functions of z, which can be taken to be increasing. 
This means that in the space of values of the SM couplings G^, variations at 
previous epochs should lie along a curve passing through Gkfl, at distances 
which increase smoothly with increasing z. 

Beyond an overall linear approximation where the (3^ are constants (for 
example, evaluated at the present value of the scalar ipo) we can consider several 
epochs through their redshifts 2;„, the corresponding values of the scalar (p{zn), 
the coupling function Bx{f{zn)) and the variations AlnGk^n- In a "local linear 
approximation" we define the coefficients Pk,n by 

AlnGk,n = dkln = Pk,nAip{Zn) ^ f5k{Zn)Aip{Zn). (9) 

We keep our assumption of proportionality of all A In Bj^ , such that A In (-^^n) — 
dkl{zn) = dkln with constant dk- In terms of the coupling functions Bk the evo- 
lution factors are defined as 

In = A\nax,n = \nBx{ip{zn)) - lnBx{(po) = PxAip{zn) (10) 

for the scenarios 2,5,6, and for the scenarios 3,4 we take 

In = Aln(</.),„ = luBni^izn)) - \^Bh{^o) = /3//A<^(z„). (11) 

3.2 Spatial variation and WEP-violating forces 

In contrast to the direct observations of time varying couplings, tests of the 
universality of free fall (weak equivalence principle, WEP) do not determine di- 
rectly the values of fundamental 'constants' or their possible variations. How- 
ever, given our basic assumption of a slow time variation, driven by a light 
scalar degree of freedom, the current limits on composition-dependent long 
range forces put bounds on the scalar couplings to different constituents of 
matter. In our language, they measure or constrain the coefficients at z = 0. 

These constraints then imply bounds on the time variation of 'constants', 
since the rate of change of the scalar is bounded above by the cosmological 
effects of its kinetic energy. If one is allowed to adjust the scalar coupling 
functions arbitrarily, the WEP force bounds on variations might be evaded, 
for example through some kind of cancellation. But in general, and for the 
unified scenarios we discuss, WEP violation places significant bounds on the 
present-day values of scalar couplings. We will follow the treatment of [131 EH] 
(see also |36j ) and consider the couplings of a scalar field to matter in the 



7 



"Einstein frame", assuming that they produce a small correction to the usual 
gravitational acceleration. 

We first briefly discuss the relation between a homogeneous cosmological 
evolution of the scalar at large scales and local variations at small scales. The 
equation of motion of the scalar <I>(t,x) = Mp(p(t,x) in a FRW cosmological 
background reads 

M + 3ff6 = -V'{<^>) - '^''^^^\ ^it,x), (12) 

mq, 

where H is the Hubble rate a/ a and ^ denotes some species of matter with 
density p^j and mass mii,($). The homogeneous solution $(t) is given by solving 

I + 3H$ = -V'i^) - '"'^^^^ pvi,(t), (13) 

where /9>i,(t, x) = p^{t) + 5p^^{t^x) and ps^ is the averaged matter density. On 
subtracting this from the general equation of motion, expanding about $ and 
making the approximations that V"{^) is negligibly small (i.e. the scalar is 
very light), that the cosmological evolution is very slow, and that the coupling 
m'^/m-i^ is much smaller than 1/Mp, we recover a Poisson-like equation 

V2(cl>_$):^!!^5p^, (14) 

applicable to small regions of space (for example the Solar System) over short 
periods of time. A more sophisticated treatment [HI [28l [29] yields similar 
results in that a local variation sourced by inhomogeneities of matter is super- 
imposed on the slow homogeneous cosmological variation. See also |30] con- 
cerning bounds on spatial variation in theories with special choices of the scalar 
couplings. 

In consequence, Newton's law for the gravitational attraction between two ^ 
particles gets modified by an additional piece which involves the dimensionless 
cosmon charge Q^i [13j 

VM = -^{l + a^). (15) 

with 

avp = 2 — = 2 (16) 

m| \ dip J 

g^ = Mp-i^. (17) 

Typically, the cosmon charge Q for two different particles is not proportional to 
their mass. For example, the proton and neutron cosmon charges Qp and Qn 
will in general not obey Qp/rup = Qn/mn, and similarly for atoms of different 
composition. As a result, two test bodies of different composition with masses 



where 
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Mb, Mc will undergo a differential acceleration towards any given gravitational 
source 

V = . (18) 

ab + ac 

We now consider how this experimentally observable quantity can be related to 
the parameters of each unified scenario, and to cosmological evolution. 

4 Apparent violation of the Weak Equivalence 
Principle 

4.1 Cosmon charge 

We may define the coupling of a dimensionless scalar ip = <I>/Mp with kinetic 
term Mp{dip)'^ /2 to a source or test body of mass Mb as 

A^.^ln^ = ^. (19) 
dip Mp Qb ^ ' 

Note that is the analogue of Pk, Eq. ([7]), if we take dimensionless mass ra- 
tios involving the Planck mass (which is fixed in the Einstein frame). Given a 
gravitational source of mass Mg separated by a distance r from two test bodies 
Mb, Mc, the static configuration for a scalar of mass m^p <^ 1/r is given by 
solving the Poisson equation with a source term proportional to A* supported 
within the source body. The resulting gradient of ip at large distance is pro- 
portional to the gravitational field of the source g, thus dp/dr = 2X'^g up to a 
sign and IS.p = /SU for small variations. H Spatial variations of couplings 
are given by Eq. ([6]) as AlnGj = difix^^, thus the "coupling coefficients" 
ki = AlnGj/Af/ used in [37] are given by ki = 2dil3x^^ ■ As discussed there, 
the resulting spatial variations within the Solar System are in principle directly 
measurable via atomic clocks, however considerable improvements in sensitiv- 
ity would be required for such measurements to compete with current WEP 
bounds. 

A test body will experience accelerations towards the source, due to the 
dependence of its mass on p, of magnitude 2g\^X'^. Thus the Eotvos parameter 
T] measured by WEP experiments is given by \13\ |35] 

??^-^ = 2^^^~2A^(A''-A^), (20) 

Qb + ttc 

where Ub^c are the accelerations towards the source of the two test masses. This 
may also be derived via "sensitivity coefficients" 

^b.c ^ dHMb/Mc) 
* dhiGi 

defined in [37], such that if^'' = Y,i ki\\^''. We then expand El 

^The factor 2 arises from the normalization of (p via the reduced Planck mass I/VSttG. 
^We use the Einstein summation convention from now on unless stated otherwise. 
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The relevant couplings A or cosmon charges Qh may be derived from the de- 
pendence of the mass of atoms on ip. It is convenient to consider the mass per 
nucleon, thus 



— In 

dip 



2/ A 



a ^ fn A '^^ 1 '^^ , D o 1 "le B B/A 

~ du, In — — -[Rz - - ] Ou, In \-Rz In d,. In , 

Mp \ 2 J rriN rriN nriN rriN Auin uin 

where uin = {nip + mn)/'2^ 5n = rUn — mp, Rz = Z/A is the proton fraction 
and B is the nuclear binding energy. 

The value of A* is typically dominated by the dependence of uin/Mp on ip. 
For unified scenarios it can be related to the (/j-dependence of I using results 
from [Ij as 

d , rriM dl , dl / dv , , \ Px , ^ 

^ = l^^^T7^-ir = dNgTT ~ { dM + 0.58— +0.35dH + 0.37ds] 23 
dl Mp dip dip \ ax J «x 

where, if ax varies with ip, we take I = In ax and dx = 1- If ax does not vary 
we replace Px/dx by Pn/dH and choose dn = 1.0 

For (a'' — A^) the leading term dln{m]y /Mp)/dip cancels. The differential 
coupling then results from the isospin- violating couplings to ^at and rUe, and 
differences in binding energy per nucleon. This last term B/A contains an 
electromagnetic contribution proportional to aZ{Z — 1)/A^/^, and a dominant 
contribution of strong nuclear forces, which to first approximation does not 
depend on A. However, in [37| the effects of subleading strong interaction 
terms in B/A, principally the surface term —asA~^^^, were included and found 
to be significant if the light quark mass m/Ac, and thus the pion mass, varies. 

To evaluate A** — A*^ we expand as Eq. ()22p using the "nuclear parameters" 
Gi = {En/niN, a, anuc/rriN}, where En = 5n - rUe, and a„„c oc ay oc is a 
parameter giving the mass scale of strong nuclear binding energy. In terms of 
these parameters we find 

A?-^ = -{R'z - Rz)-Q^ (^) - - ■ (24) 

With ARz = R^z-^z this yields 

A^7 = -ARz^, (25) 

while the dependence on the other nuclear parameters is given by the binding 
energy and is discussed in [37] . 



*If we take a alone to vary, the only contribution to A** is from the electromagnetic contribu- 
tions to the proton mass and nuclear binding energies, hence we evaluate this case separately. 
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In a given unified scenario the dependences of Gi on {p variation arise via 

d\nGi dl 

and one finds 

1 29 

dEn ^ —dx + IAUh - 0.82ds, 

ax 

, 80a , 43 a , 257 a , 
^"=27^^^ + 27 2^^^ + ^2;^^^' 

d.,nu. - -O.g^MlVAe) ^ 063,^ _ , ^ , (27) 
A In ax ax 

where Eqs. (24), (25) and (44) of [Ij were used for the dependence of En, anuc, 
and rriN/M-p on fundamental parameters. 

In any given unified model the couplings A**, A*''^ are proportional to the 
coupling of if to the fundamental parameter. For a nonzero j3x {(3h) we thus 
have r] oc {Ph) fundamental coupling at the present value of (p is 

directly bounded. 

4.2 Unification factor 

We would like to compare the differential acceleration ij with the time variation 
of couplings. In order to arrive at the simple formula ([T|), we first consider a 
"fiducial unified scenario" for which the unification factor satisfies F = ARz 
in this equation. Other unified scenarios can then be characterized by the 
difference of F/ARz from unity. 

As fiducial unified scenario we consider a variation of ax in a non-SUSY 
GUT (ax — 1/40), accompanied by a variation of {(p) such that the ratio 
between the Fermi scale {4>) and the QCD scale Ag remains constant. This 
corresponds to the unification coefficients 

2tt 

du = 0, dx = 1, dn = ^ — ^ 36, ds = 0. (28) 

Tax 

We also neglect for the fiducial calculation the contributions from nuclear bind- 
ing energy. These can be incorporated into the unification factor via 



F = Cgut (^ARz + DgutA^^^^ + D'auT^A-^'""^ 



(29) 



where Cgut accounts for the source coupling and value oi dEn, and the terms 
with prefactors Dgut and Dq^j, account for nuclear Coulomb self-energy and 
nuclear surface energy, respectively. 

For the fiducial calculation the preceding relations yield 



b-c o^sA.. -En dl „0.90^„ -En, ^r.^, , fdlnax 



2X'ARz '-^dE„^ = 2^ARz ^(-0.97d,) , (30) 

niN oip ax mx \ 0(p J 



2 
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using the relation dEn — —0-97da that holds if fermion masses do not vary 
relative to Ac. Now writing 

9 In ax 1 dlnax d~^ f 2877q!x 1 /d\ 

d(p 1^ dt if) \a J 67 \a J ' 

and using the cosmological relation 99^ = SH^Qhil- + Wh) we arrive at 

Here fi/j denotes the fraction of the energy density of the scalar field as compared 
to the total energy density of the Universe. The equation of state Wh = {T — 
V)/{T + V) can range between -1 and +1 for a single scalar field. If the scalar 
is responsible for the dark energy in the Universe, one needs Q.h ~ 0.73 and 
1 + ^ 0.1 [38] . The relation (132]) is also independent of the value of ax, as 
long as Eq. ([28]) is satisfied. 

We next consider the contributions from Coulomb self-energy and nuclear 
surface energy. The ratio of the Coulomb term to the Ai?^ term in 77 yields in 
general 

, 80a A I 43 a J , 257 « j 

DauT = -^4^ - -0.91 ^^^ + ^^'^^ + ^^^^ , (33) 
-En dEn -1.29a^^dx + 1.41^^ - 0.82^5 

where ac — 0.71 MeV. In our fiducial unified scenario the strong binding energy 
coefficient anuc/fnN does not vary, thus Dqut — o,c/{0.97En) — 0.94 and 
E^GUT ~ 0- More generally we have 



j^, _ as dg^nuc ^ g O.eSax^dx - 0.7dH + 0.4dg 

-En dEn ~ ' -1.29a]^^dx + lAldH -0.82ds' 



where as — 17.8 MeV. 

By definition the fiducial scenario has Cqut = 1- In a general unified 
scenario we insert the appropriate values for the unification coefficients dk, 

CauT= '^^ > f^)', (35) 

where A^q^, d_E„(o) ^^a(o) the values for the fiducial scenario. In Tabled] 
we display CcuT, Dqut and D'qjjj, for the various unification scenarios. We 
also display the values of E obtained by combining the different ingredients in 
Eq. ()29p . and considering Be and Ti test masses (ARz — —0.015, AZ(Z — 
1)^-4/3 ^ _2,o, ^ Q 21). The very large value of E for Scenario 3 

reflects an accidental cancellation of contributions to dlna/dl which we do not 
consider to be typical. 

4.3 WEP limits on variations at the present epoch 

The time evolution of a scalar contributes to the effective cosmological equation 
of state via its kinetic term to bound at present, in terms of the energy density 
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Scenario 


' 7ax 


2 


3 


4 


5, 7 = 42 


6, 7 = 70 


6, 7 = 25 


Cgut 


1 


43 


-3900 


-69 


-8.8 


-9.0 


19 


Dgut 


0.94 


0.015 


-0.001 


-0.020 


-0.11 


-0.12 


0.048 


n' 





11.1 


11.3 


11.5 


12.7 


12.8 


10.7 


F (Be-Ti) 


-1.90 


95 


-9000 


-165 


-25 


-26 


41 



Table 1: Values of Cgut, ^gut "fiducial scenario" and for unified 

scenarios considered in [1]; we also give values of F for a WEP experiment 
using Be-Ti masses. 



fraction fi/j and equation of state Wh of quintessence. Previously [35] we took 
a conservative upper limit of ^^^(l -|- Wh) < 0.16 leading to a "speed limit" 
ip < ifmax — 5 X 10~^^ for a Hubble rate of 7 x 10~^^ y~^- Considering recent 
improvements in cosmological measurements we may now take Wh < —0.9 and 
— 0.73 at the present epoch (if the evolution of the field has not accelerated 
at late times), thus ^/Hq < 0.47 and we have 

Vimax =^3.5 X 10-"y-i. (36) 

Combining this with the upper limit from WEP on the fundamental coupling 
f3x or Ph in each unified model, we bound the present rate of variation of the 
fundamental parameter ax or (i;^) /Mx- The bounds can be extrapolated to past 
epochs by expanding the coupling functions and the scalar time evolution about 
the present values, if these are smooth and not too rapidly-varying functions. 

The experiment setting the tightest limits on scalar couplings has the 
result 

rj = (0.3 ± 1.8) X 10-^3 (37) 

for test bodies of Be and Ti composition, where the gravitational source is taken 
to be the Earth. 



Variation of a alone Here we parameterize the variation as A In a = PaA(p, 
and we have for the source coupling 

= f 3 X 10-4 + 8 X W-^Rz - 1/2) + — ^^^^l Pa (38) 

where the coefficients arise from the contributions of electromagnetic self-energy 
to mjv, 6iy and B/A respectively. We have dEn — —0.97, and clearly also da = 1, 
da,nuc = 0. The bound arising from Eq. ([37]) . treating the Earth as composed of 
iron or silicon (with similar results), is /3a < 1.8 x 10"^ at la. Thus the present 
rate of variation is limited to \a/a\ < 6.2 x 10"^^ y"^; a somewhat looser bound 
than from atomic clocks. 



WEP bounds for unified models If the unified coupling ax varies then 
our fundamental coupling is [3x- For Scenario 2 we find the bound 

^1 < 6.2 X lO^^^Q^ (Icr) 
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thus the present rate of variation is hmited to 

\ax/oix\ < 2.8 X W'^^axy'^, (39) 

equalhng 1.2 x IQ-'^'^ y-^ for SUSY theories {ax ^ 1/24) or 6.9 x lO-^'^y-^ 
for non-SUSY {ax — 1/40). A major contribution to rj is the strong nuclear 
binding energy with coefficient as/vriM, which is sensitive to variation of the 
Ught quark mass relative to the QCD scale. 

For Scenarios 3 and 4 we take instead a fundamental coupling (3h describing 
the variation of {(())/ Mx (and fh/Mx in the SUSY case). In Scenario 3 the 
magnitude of I3h is limited to 9.5 x 10^^ (la) and time variation is bounded 
by |9f ln(((^)/Mx)| < 3.3 x 10~^^. For Scenario 4 the bounds are similar with 
Ps < 1-0 X 10-5 and \dtln{{4>)/Mx)\ < 3.5 x lO'^^. 

For Scenarios 5 and 6 we take again a unified coupling Px- in Scenario 5 
with 7 = 42 we find an upper limit \f3x\ < 3.6 x 10"^, thus axjoix is bounded 
at 1.3 X lO^i'^y^i at \a. 

In Scenario 6 with 7 = 70 the relevant \o bounds are still tighter with 
1/9x1 < 2.3 X lO"'^ and |dx/ax| < 8.1 x lO'^^y'^ With 7 = 25 here we obtain 
1/3x1 < 3.1 X 10"'^ and |dx/ax| < 1.1 x lO^^'^y-^ 

Cancellations in WEP bounds For particular values of 7 in Scenarios 5 
and 6 there may be cancellation between different contributions to r/, thus A;, 
and Ac may happen to be nearly equal for the test bodies considered. Without 
SUSY (Scenario 5) this happens for 7 ~ 35, approximately the value 807r/7 at 
which the variation of ((/>)/Ac also vanishes. With SUSY (Scenario 6) and dg 
set equal to du^ the analogous value of 7 for which A In (0)/ Ac = is 167r ~ 50. 
Near these values, contributions from variation of fermion masses relative to Ac 
become small, leaving the intrinsically smaller effect of varying a. In principle 
one might fine-tune 7 to produce a yet more dramatic cancellation in any given 
experiment. However the cancellation cannot be exact for all pairs of test 
bodies, thus bounds on /3x from all WEP experiments cannot be removed by a 
fine-tuned choice of 7, only slightly weakened. 

4.4 Bounds on present-day variation and testing unified 
scenarios 

Within our theoretical framework there exist three distinct ways to bound or 
measure the present-day rate of variation of fundamental parameters. The 
first is a direct measurement, for instance atomic clock experiments. If one or 
more nonzero variations is found in this way, bounds on unified models may 
immediately be set. The second method is by combining information on the 
size of scalar field couplings from WEP tests (Section 13. 2p with a cosmological 
upper bound on the kinetic energy of scalar fields [13l [35]. Such bounds on 
scalar couplings will depend on the choice of unified model and in general will 
be independent of those derived from atomic clocks. 

Thirdly, under the assumption of a monotonic variation (that also does not 
significantly accelerate with time), we may convert any "historic" bound on the 
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net variation of a fundamental parameter since a given time, into a bound on 
the present rate of variation: 

\Gk\<{to-tn)-'\Gk{to)-Gk{tn)\=^-^^, tn<to, (40) 

where Iq denotes the present time. 

For any given unified model of time variations, the three methods to bound 
present-day evolution will have different sensitivity. Therefore if one method 
gives a clear nonzero variation we would (in some cases) be able to distinguish 
between models, given that the other bounds are still consistent with zero. 
To give a simple example, the direct detection of a nonzero time variation 
in atomic clocks near the present upper bound would immediately rule out a 
large class of models that cannot account for such a variation without leading 
to WEP violation above current limits, and would also rule out models in 
which such nonzero variations extrapolated to past epochs t„ < to would exceed 
observational bounds. 

However, such inferences do not function equally in all directions. A nonzero 
finding of WEP-violating differential acceleration would indicate nontrivial scalar 
couplings, but need not imply nonzero time variation since the rate of change 
of the scalar is not bounded below. Also, a nonzero variation at some past 
epoch tn would not necessarily imply a lower bound to the present-day rate 
of variation or size of scalar couplings, since the variation could have slowed 
substantially since then (either due to nonlinear scalar evolution or a nonlinear 
coupling function). With the assumption of a smooth and monotonic varia- 
tion of the scalar field and its coupling functions, one could estimate, for any 
given unified model, where the first signals of present-day or recent variation 
are expected to appear. 

At present these methods give null results up to redshifts about 0.8, but if a 
nonzero time variation exists, we can determine for each unified scenario which 
observational method is most sensitive. Thus if a nonzero signal of recent vari- 
ation arises it may be used to distinguish between models. We assume for this 
purpose an approximately linear variation over recent cosmological times, thus 
measurements of absolute variation at nonzero redshift z imply time derivatives 

dlnX _ AlnX{z) 
dt to — t{z) 

Here X is the fundamental varying parameter: we consider first X = a, the 
only varying parameter is a; in scenario 1, X = G^m\, in scenarios 2, 5 and 
6, X = ax and in scenarios 3 and 4, X = {(j))/Mx- Then Table [2] gives the 
precision of bounds on time derivatives for the unified scenarios we consider. H 
As in [1] we take the Oklo bound as applying directly to the variation of a, and 
increase its uncertainty by a factor 3 to account for possible cancellations when 
other parameters also vary. We present the Rosenband et al. [40j Al/Hg ion 
clock bound separately to illustrate to what extent it improves over previous 
atomic clock results. 

^Except scenario 1 (varying Gn) which is probed by a quite different set of measurements. 
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(41) 



Scenario 


X 


Clocks 


Al/Hg 


WEP 


Oklo a 


Meteorite 


Astro 


a only 


a 


0.13 (a) 


0.023 


6.2 


0.033 


0.32 


0.44 {y) 


2 


ax 


0.074 (^) 


0.027 


0.007 


0.12 


0.015 


0.006 (NH3) 


2S 


ctx 


0.12 (/i) 


0.044 


0.012 


0.19 


0.026 


0.010 (NH3) 


3 


{ct>)/Mx 


2.6 (m) 


12.4 


0.33 


54 


0.53 


0.22 (NH3) 


4 


{(p)/Mx 


6.2 (^) 


1.78 


0.35 


7.7 


1.2 


0.51 (NH3) 


5, 7 = 42 


ax 


0.32 (a) 


0.024 


0.013 


0.11 


0.069 


0.035 (NH3) 


6, 7 = 70 


ax 


0.21 (a) 


0.016 


0.008 


0.070 


0.049 


0.025 (NH3) 


6, 7 = 25 


ax 


0.25 (^) 


0.027 


0.011 


0.12 


0.056 


0.021 (NH3) 



Table 2: Competing bounds on recent {z < 0.8) time variations in unified 
scenarios. For each scenario we give la uncertainties of bounds on d{\nX)/dt 
in units 10~^^y~^, where X is the appropriate fundamental parameter. The 
Oklo bound is rescaled in Scenarios 2-6 as explained in the main text. The 
column "Clocks" indicates whether a oi n gives the stronger bound; the recent 
Al/Hg limit [40] is given a separate column. The column "Astro" indicates 
which measurements of astrophysical spectra are currently most sensitive in 
each scenario. 

Extending the methods of this section beyond z ~ 0.5 becomes questionable. 
One could use linearity in ln(l + z) instead of t, but even this improvement may 
lead to unreliable extrapolations for models with a particular dynamics of the 
scalar field, such as crossover quintessence. An alternative approach to relating 
present-day variation to cosmological history in previous epochs, under certain 
assumptions on the scalar evolution, is given in [41j . 

5 Monotonic cosmon evolution and crossover 
quintessence 

5.1 Crossover quintessence 

We now consider two illustrative models which predict the evolution of the 
quintessence field ip{t), consistently with other cosmological bounds, given a 
small number of adjustable parameters. Our aim is to see to what extent such 
models can be consistent with the behaviour of time variations that we have 
outlined. 

Our first class of models is "crossover quintessence" [42^ l43l [5]: here the 
scalar field follows tracking solutions [3l SI] at large redshift. In this early 
epoch the equation of state Wh = {T — V)/{T + V), where the kinetic energy 
T = Mp93^/2 and potential energy is V{ip), is equal to that of the dominant 
energy component (matter or radiation). At some intermediate redshift before 
the onset of acceleration, the time evolution of the cosmon slows down. In 
consequence, there is a crossover to a negative equation of state and the fraction 
of energy density due to the scalar begins to grow. In recent epochs the field 
has an effective equation of state Wh ^ —1- We do not aim to build and solve 
detailed models of this type, but rather estimate general properties of the scalar 
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Scenario 


h/k 


h/h 


h/h 







0.00 ±0.01 


0.13 ±0.13 


0.31 ±0.18 


-600 ± 4000 


2 


-0.10 ±0.12 


0.04 ± 0.04 


-1.59 ±0.88 


-380 ± 230 


3 


-0.11 ±0.13 


0.04 ±0.03 


-0.12 ±0.91 


-390 ± 240 


4 


-0.04 ±0.10 


0.04 ±0.03 


0.97 ±0.64 


-380 ± 260 


5, 7 = 42 


0.00 ±0.03 


0.04 ±0.03 


0.41 ±0.19 


-280 ± 191 


without BBN 


11 






-66 ± 165 


6, 7 = 70 


0.00 ±0.03 


0.02 ± 0.02 


0.39 ±0.17 


-275 ± 150 


without BBN 








-69 ± 139 


6, 7 = 25 


-0.04 ± 0.06 


0.02 ± 0.04 


-1.04 ±0.49 


-343 ± 142 


with Wendt 


0.03 ±0.05 


-0.01 ±0.03 


0.70 ±0.52 


231 ± 163 



Table 3: Ratios of the evolution factors from observations. 



evolution. 

We begin with a caricature of a scaling quintessence or early dark energy 
model, in which the dark energy equation of state is constant at late times at 
Who- Above some given redshift z+ the equation of state crosses over to the 
scaling condition lu/j = in the matter dominated era; then for z > z^q, before 
matter-radiation equality, we again have scaling through Wh = 1/3. Then the 
general relation (a = (1 ± z)^^) 

^ = -3(1 +w{a)) (42) 
a ma 

may be used to find the matter, radiation and dark energy densities over cosmo- 
logical time. We estimate the scalar kinetic energy via M'^tp^ /2 = ph{^ + Wh)/'2. 
and can thus integrate dip /da = ip/aH from the present back to any previous 
redshift. The initial conditions are set by specifying the present densities of 
matter, radiation and dark energy and the model parameters Who and 

Assuming a constant coupling 5 to the fundamental varying parameter, 
usually ax, the variation is given by 

^\nax{z) = 5{p{z)-p{{})). (43) 

This ansatz implies a monotonic evolution, since the increase of p with time is 
monotonic for crossover quintessence. We discussed the viability of monotonic 
evolution in [T] where a first judgement can be made by inspection of figs 1-6 
for the various unification scenarios, or Table [6] in the appendix of the present 
paper. 

To allow us to compare easily with the observational results, we observe 
that in Table [6] the constraints for h are considerably weaker than those for 
the remaining evolution factors. Furthermore, the one-sigma range for Z4 is 
nonzero for all scenarios. Hence we shall compare observational and theoretical 
values for the ratios /1//4, h/li-, h/h and h/h- We note the opposite sign of ^3 
and I4 for scenarios 2 and 3, and scenario 6 with 7 = 25, which disfavours any 
monotonic evolution for these scenarios. The averaged observational values in 
each epoch are given in Table [3l 



17 



Who 


z+ 


h/k 


k/k 


h/k 




-0.95 


3 


0.12 


0.38 


0.67 


26 


-0.99 


3 


0.09 


0.28 


0.55 


38 


-0.9999 


3 


0.02 


0.10 


0.37 


54 


-0.95 


7 


0.15 


0.47 


0.80 


13 


-0.99 


7 


0.15 


0.47 


0.80 


24 


-0.9999 


7 


0.08 


0.27 


0.52 


120 


-0.999999 


7 


0.02 


0.09 


0.34 


180 



Table 4: Ratios of evolution factors from crossover quintessence. The li are 
evaluated by averaging over the variations evaluated at the same redshift as 
the data in each epochs (weighting by the number of absorption systems if 
appropriate) . 

Note that a constant coupling 5 drops out of the ratios li/lj. Thus we may 
probe cosmon/quintessence evolution directly without knowing the absolute size 
of the coupling, only assuming its (approximate) constancy over the relevant 
range of evolution. Due to its monotonic evolution, crossover quintessence 
cannot give negative ratios U/lj. Hence it cannot be a good fit to the ^Li 
abundance within the unification scenarios 2 to 6 that we consider. This reflects 
the tension between the ^Li abundance and a positive variation of fi discussed 
in Section 4.4 of [l]. 

In Table m we display the values of the ratios h/l^ and li/l^ expected from 
crossover quintessence for various values of the parameters Who and . Consid- 
ering the ^Li abundance to be affected by astrophysical systematics in scenarios 

5 and 6 (7 = 70), or using the null result of [27j for Aln/i at intermediate 
redshift for scenario 6 (7 = 25 "with Wendt"), we conclude that crossover 
quintessence could, in pinciple, reconcile a coupling variation of the claimed 
size in epochs 3 and 4 with the bounds from late cosmology, i.e. epochs 1 and 
2. This is due to the "slowing down" of the cosmon evolution, as noted in |14j . 
Values of the present equation of state Who quite close to —1 would be required, 
however. An observation of coupling variations would put strong bounds on the 
dynamics of the cosmon field and provide an independent source of information 
about the properties of dark energy. 

Note that observational probes of dark energy would not give results for 
Who that coincide with the values that we take in our model. Such probes do 
not actually measure the present-day equation of state, rather they estimate wq 
by extrapolating from past epochs under some parameterization. 

6 Models with growing neutrinos and oscillating 
variation 

In this section we investigate models that do not obey the proportionality of all 
coupling variations for all redshifts and do not show a monotonic evolution of 
the cosmon field. A systematic analysis of all such models seems difficult, and 
we concentrate on a specific example. 
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Growing neutrino models \2Q\ [2T] explain the value of today's dark energy 
density by the "principle of cosmological selection". The present fraction of 
dark energy, Q^, is set by a dynamical mechanism. As soon as the neutrinos 
become non-relativistic, their coupling to the cosmon triggers an effective stop 
(or substantial slowing) of the evolution of the cosmon. Before this event, the 
quintessence field rolls down an exponential potential and thus follows the track- 
ing behaviour described in the preceding section. This mechanism requires a 
neutrino mass that depends on the cosmon field and grows in the course of the 
cosmological evolution. The present dark energy density, p^Q, can be expressed 
in terms of the average present neutrino mass, m,y(to), and a dimensionless 
parameter 7 of order unity [20], 

ip,o)'/^ = l.07{l^y\o~^eY. (44) 

We follow again our simple proportionality assumption, namely that the 
cosmon evolution produces a variation Alnax{z) = 5{ip{z) — (^(0)), with a 
proportional variation for other couplings according to the unification scenario. 
This is the only contribution to the variation of the unified coupling ax and 
Mx/Mp. However, a new ingredient is an additional variation of the Higgs 
v.e.v. (i;^) with respect to the Planck mass, which only becomes relevant at late 
time [21]. It is due to the effect of a changing weak triplet operator on the 
v.e.v. of the Higgs doublet. If the dominant contribution to the neutrino mass 
arises from the "cascade mechanism" (or "induced triplet mechanism" ) via the 
expectation value of this triplet, this changing triplet value is directly related 
to the growing neutrino mass |2lj. 



6.1 Stopping and scaling growing neutrinos 

We consider two models, with slightly different functional dependence of the 
Higgs v.e.v. and neutrino mass on the scalar field. In the first, where the cos- 
mon asymptotically approaches a constant value ("stopping growing neutrino 
model") [2IJ, the additional Higgs variation is given according to 

^'^^(z) = ^(l-/^(z))-°•^ (45) 



Mx 
where 

1 - exp{-e{ip{z) - ipt)) 

Here, (pt ~ 27.6 is the asymptotic value (choosing the parameter a = 10 in 
the exponential potential [21j). For illustration we take the set of parameters 
given in [21], e = —0.05, while H is set by demanding the Higgs v.e.v. being 
consistent with measurements today, {(j)){z = 0) = 175 GeV. We set Rq = 10"'', 
however in general we only require R{z = 0) ^ 1. The resulting variations are 
shown in Fig. [TJ 

The combined variation is different for each of the unified scenarios. For all 
unified parameters except (</>) /Mx we still have a proportionality for the vari- 
ations at all z, since the variations are proportional to ip. The variations due 
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0.5 1 2 4 6 9 



A'ln<(^>/Mx[10-4] 

0.5 



1 



2 



4 



6 



9 



z 



-0.05 



-0.10 



-0.15 



-0.20 



Figure 1: Variations in the stopping growing neutrino model of [21j. Left: 
evolution of the cosmon field. Right: additional variation of the Higgs v.e.v. 
according to Eq. (jl5|) . 

to the direct coupling (jl3]) can be described by our method of evolution factors 
In, even though the In need not be strictly monotonic due to the oscillations 
in ip. However, for {(p)/Mx we now have one variation linearly proportional 
to the variation of (p, Eq. (143p and an additional one with a nonlinear depen- 
dence, Eq. (j45ti46p . A simple treatment with common evolution factors for all 
variations will no longer be applicable. Due to the additional (/j-dependence of 
{4>)/Mx we may have separate evolution factors for {(f)) /Mx, different from the 
(common) evolution factors for the other couplings. 

For example, the "linear contribution" (I43p may dominate at BBN and in- 
duce a positive Iq common for all couplings. In the range z < 10 the "non-linear 
contribution" (j45p could be more important, leading to effectively negative ^3^4 
for {(f)) /Mx- (Such an effect could, in principle, relieve the tension between ^Li 
and a positive ^u-variation at high z.) In practice, we calculate Alna^ and 
A In at each epoch directly from the model, and extract the varying cou- 
plings and observables as explained in Section 3 of P]. Then we may search for 
a set of parameters 5, Rq which minimizes the for all measured variations. 

The stopping growing neutrino model has an oscillation in (i;^) that grows 
both in frequency and amplitude at late times as approaches its asymptotic 
value. Such oscillations must not be too strong as measurements between z = 2 
and today would measure a high rate of change. The oscillation may be made 
arbitrarily small by choosing small Rq. The restrictions from the low-z epochs 
are actually so strong that to a good approximation the non-linear contribution 
~ Rq can be neglected. However, the linear variation (I43p is independent of 
Rq. It can be described by our method of evolution factors and yields for our 
set of parameters the ratios 



h/U = 0.008, 

h/h = 0.09, 

h/U = 0.44, 

h/U = 175. 



(47) 
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Figure 2: Variations in the scaling growing neutrino model of [20]. Left: evolu- 
tion of the cosmon field. Right: additional variation of the Higgs v.e.v. accord- 
ing to Eq. dSD. 

Comparing this with the numbers given in Table [3] shows that this model nat- 
urally yields evolution factors which are of the correct order of magnitude. We 
emphasize that no new parameter has been introduced for this purpose. 

The second growing neutrino model [20] does not lead to an asymptotically 
constant ip. Now the coupling of the neutrino to the cosmon ip is given by a 
constant /3, according to 

rriu = rh^e-^'^. (48) 

This "scaling growing neutrino model" leads in the future to a scaling solution 
with a constant ratio between the neutrino and cosmon contributions to the 
energy density. 

With the choice of parameters /3 = —52, a = 10 and nij^fi = 2.3 eV, and 
given the triplet mechanism of [21], the Higgs v.e.v. varies as Eq. (j45p . where 
now R is given by 

R{z) = Roe-^^^'\ (49) 

Here the Higgs oscillations remain comparatively small in amplitude, while the 
absolute value of ((/>) grows overall with time: see Fig. [2] with the parameter 
choice Rq = 10^^. Most of the additional variation of the Higgs v.e.v. occurs 
at later epochs, z < 2, thus recent observational bounds rule out any signifi- 
cant additional growth in ((/>). We considered fitting the observational values 
excluding BBN, as a function of the model parameters 6 and Rq, and we find 
always that the value of Rq at the minimum of is unobservably small. 

6.2 Global fits to growing matter models 

Each growing neutrino model contains a few parameters that determine the 
cosmological evolution of the cosmon ip, and the Higgs v.e.v. ((/>). Two coupling 
parameters give respectively the relative strength of variation of ax with p, and 
the relative strength of the additional variation of the Higgs v.e.v. due to the 
varying triplet. For each example of cosmological evolution we may calculate 
the observables directly in terms of the two coupling parameters and make a 
global fit for their values. In performing the fit we take the 125 systems of the 
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Murphy et al. a determination [22] within Epochs 3 and 4 (Eq. (9) of T]) and 
further spht them into 5 subsamples each with 25 absorption systems, since the 
data set extends over a wide range of redshift where there may be significant 
oscillations. 

For the global fits, we take the scaling growing neutrino model [2D] with 
P = —52, a = 10, m^fi = 2.3 eV. With zero variation at all times (no degrees 
of freedom), we find = 3.2 including "^Li at BBN, and = 2.4 neglecting 
^Li. The results of the best fits with varying couplings are given in Table [5l No 



Scenario 


5 X 10^ 


Ro 


x' 


Ax' 


2 


-0.019 


0.045 


3.1 


0.2 


2 (no ■^Li) 


-0.040 





2.1 


0.3 


3 


1.6 





2.8 


0.4 


3 (no ■^Li) 


1.6 





2.0 


0.4 


4 


3.8 





2.8 


0.4 


4 (no ■^Li) 


3.8 





2.0 


0.4 


5.42 


0.30 





2.0 


1.3 


5.42 (no ■^Li) 


0.24 





1.9 


0.5 


6.70 


0.20 





1.9 


1.3 


6.70 (no ■^Li) 


0.16 





1.9 


0.5 


6.25 


0.18 


0.090 


2.7 


0.5 


6.25 (no ■^Li) 


-0.055 


0.061 


2.2 


0.2 



Table 5: Fitting parameters and minimal values for the different unification 
scenarios for best fit to the scaling growing neutrino model [20j. The last column 
gives the increase in x' produced when 6 and Rq are forced to vanish, i.e. for 
zero variation. 

convincing evidence is found for coupling variations within this model. We have 
investigated some other choices of parameters and also the stopping growing 
neutrino model, without a substantial change in the overall situation. In view 
of the unsettled status of the observational data it seems premature to make a 
systematic scan in parameter space. Our invesigation demonstrates, however, 
how a clear positive signal for a coupling variation could restrict the parameter 
space for quintessence models. 

7 Conclusions 

This paper demonstrates how a clear observation of time variation of funda- 
mental couplings would not only rule out a constant dark energy, but also put 
important constraints on the time evolution of a dynamical dark energy or 
quintessence. 

In Sec. m we have seen how the comparison of a varying coupling with the 
bounds from tests of the equivalence principle can put a lower nonzero bound 
on the combination 0,^(1 + 11)^), according to 

^hil + Wh) > 3.8 X 10i8F(at lnQ)2r/-L , (50) 
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with dt In o. in units of y , ctnci where ^max — 

1.8 X 10~^3 is the current ex- 
perimental hmit on the differential acceleration of two test bodies of different 
composition. Thus, if |f?tlna| is nonzero and not too small, Wh cannot be ar- 
bitrarily close to —1 (a cosmological constant); nor can the contribution of the 
scalar to the dark energy density be insignificant. 

The precise bound depends on the "unification factor" F, which differs be- 
tween different scenarios of unification and also depends on the composition of 
the experimental test bodies. We find 1 < -F < few x 10^ for the Be-Ti masses 
used for the best current bound on r/ and for representative cases of unified 
scenario. 

Conversely, assuming that the scalar field responsible for the varying cou- 
plings is the cosmon, whose potential and kinetic energy account for a dynam- 
ical dark energy in the Universe, the bounds on ij can be used to set bounds 
on the time variation of couplings in the present epoch. For the different uni- 
fied scenarios we compare the relative sensitivity of these bounds as compared 
to laboratory measurements or bounds from the Oklo natural reactor and the 
composition of meteorites in Table [2j 

For a given unified scenario, the bounds on the time variation of various 
couplings in different cosmological epochs strongly restrict the possible time 
evolution of the cosmon field, once at least one irrefutable observation of some 
coupling variation at some redshift becomes available. We have demonstrated 
this by an analysis that implicitly assumes a nonzero variation, considering both 
general features and specific quintessence models. We are aware that the actual 
values for the evolution factors l{z) from this analysis may be premature, since 
the observational situation is unclear and on moving grounds. For example, 
taking the recent reanalysis of the variation of the proton to electron mass ra- 
tio fj. in Ref. [l5] instead of the results in Ref. [23] used in this paper, would 
strongly influence the values of the evolution factors. We have demonstrated 
this in a somewhat different way by investigating the change in the evolution 
factors if some claimed observations of varying couplings are omitted. Needless 
to say that without a clear signal any analysis may become obsolete. For the 
time being our analysis remains a useful tool for the comparison of the rela- 
tive sensitivity of different experiments and observations, and for a judgement 
about mutual consistency of different claimed variations and bounds from other 
observations. 

Acknowledgements 

We acknowledge useful discussions with M. Doran, G. Robbers, M. Pospelov, 
J. Donoghue and P. Avelino, and correspondence with V. Flambaum. T. D. 
thanks the Perimeter Institute for hospitality during the workshop "Search for 
Variations in Fundamental Couplings and Mass Scales". T. D. is supported by 
the Impuls- and Vernetzungsfond der Helmholtz-Gesellschaft. 



23 



Appendix: Evolution factors 



Here we reproduce for convenience of the reader the table of evolution factors 
for the unification scenarios and epochs considered in [T]. 



Epoch 


i 


O 

Z 


6 


4 


r 



D 




0.14 


0.53 


1.6 


3.8 


10^ 


IQio 


Scenario 


h X 10** 


h X 10** 


h X 10^ 


U X 10^ 


h X lO-* 


Iq X 10^ 


a only 


-0.01 ±0.06 


-1.1 ± 1.0 


-0.26 ±0.10 


-0.85 ±0.37 


-150 ±350 


5 ±34 


2 


-0.1 ±0.1 


0.04 ± 0.03 


-0.15 ±0.08 


0.10 ±0.03 


0.9 ± 14 


-0.37 ±0.20 


3 


4.1 ±4.8 


-1.5 ± 1.2 


0.42 ±3.3 


-3.6 ±0.9 


69 ± 920 


14 ±8 


4 


3.9 ±8.5 


-3.4 ±2.7 


-8.4 ±5.1 


-8.7 ±2.1 


31 ± 450 


33 ±21 


5, 


-0.02 ±0.18 


-0.24 ±0.18 


-0.25 ±0.10 


-0.61 ±0.13 


0.6 ±8.6 


1.7± 1.1 


(7 = 42) 












[0.4 ± 1.0] 


6, 


-0.02 ±0.12 


-0.10 ±0.07 


-0.17 ±0.07 


-0.44 ±0.10 


0.3 ±5.0 


1.2 ±0.6 


(7 = 70) 












[0.3 ± 0.6] 


6, 


-0.12 ±0.18 


0.04 ±0.12 


-0.30 ±0.11 


0.29 ± 0.08 


0.7 ± 10 


-1±0.3 


(7 = 25) 








[-0.43 ± 0.28] 







Table 6: Redshifts and evolution factors for each epoch, for each scenario. In 
the first row the values of In give the fractional variation of a; in Scenarios 2, 
5 and 6 that of ax] and in 3 and 4 that of {(f)) /Mx- Values in brackets give, 
for BBN (/g), the evolution factors neglecting ''Li; or for /4, the evolution factor 
with the Reinhold values [23] substituted by the Wendt value [27j . 
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